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We have studied constraints on the equation of state, w, and speed of sound, c s , of the dark energy 
from a joint analysis of data from the cosmic microwave background, large scale structure and type- 
la supernovae. We find that current observations have no significant sensitivity to c s . However, there 
is a slight difference between models in which there are no dark energy perturbations and models in 
which dark energy behaves as a fluid. Assuming that there are no dark energy perturbations shifts 
the allowed region for w to slightly higher values. At present models with and without dark energy 
perturbations provide roughly equally good fits to observations, but the difference is potentially 
important for future parameter estimations. Finally, we have also performed error forecasts for 
future measurements of c s . 



PACS numbers: 



I. INTRODUCTION 



The discovery in 1998 from observations of type la su- 
pernovae [J, that the universal expansion is currently 
accelerating was a spectacular result. The finding has 
since been confirmed by observations of the Cosmic Mi- 
crowave Background (CMB) [101 and the large scale 
structure (LSS) of the universe 0,0,0 

One possible explanation is that the energy density of 
the universe is dominated by dark energy with a negative 
equation of state. The simplest possibility is the cosmo- 
logical constant which has P — wp with w = —1 at all 
times. However, since the cosmological constant has a 
magnitude completely different from theoretical expecta- 
tions one is naturally led to consider other explanations 
for the dark energy. 

A light scalar field rolling in a very flat potential would 
for instance have a strongly negative equation of state, 
and would in the limit of a completely flat potential lead 
to w = —1 [1, 0, El • Such models are generically known 
as quintessence models. The scalar field is usually as- 
sumed to be minimally coupled to matter, but very in- 
teresting effects can occur if this assumption is relaxed 
(see for instance prjp. 

In general such models would also require fine tuning 
in order to achieve fix ~ flm , where fix and fl m are the 
dark energy and matter densities at present. However, 
by coupling quintessence to matter and radiation it is 
possible to achieve a tracking behavior of the scalar field 
so that fix ~ flm comes out naturall y o f the evolution 
equation for the scalar field [3 [3 Q [3 13 13 13 13 

Many other possibilities have been considered, like k- 
essence, which is essentially a scalar field with a non- 
standard kinetic term [3 III 13 13 13 13 El- It is also 
possible, although not without problems, to construct 
models which have w < — 1, the so-called phantom energy 
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Finally, there are even more exotic models where the 
cosmological acceleration is not provided by dark energy, 
but rather by a modification of the Friedman equation 
due to modifications of gravity on large scales [4J, |45| . 



If dark energy is a fluid, i.e. a component with en- 
ergy density and pressure, it can in general be charac- 
terized by one parameter in addition to the equation of 
state, w = Pj p, namely the speed of sound, c s . Scalar 
field models with standard kinetic terms all have c s = 1 
which is a solid prediction. On the other hand, mod- 
els such as fc-essence can have an effective sound speed 
which at times is c s 3> 1 (note that this in fact does 
not violate causality). Models where the acceleration is 
due to a modification of gravity have no definable sound 
speed, and the "dark energy" component does not have 
perturbations. Measuring the effective sound speed of 
dark energy would therefore provide crucial information 
which is complementary to the measurement of the equa- 
tion of state. Note that another possibility is to discuss 
the dark energy parameters purely in terms of observ- 
ables instead. This approach was for instance adobted in 
the recent paper ^(j . 



In the present paper we discuss current constraints on 
the sound speed of dark energy as well as the prospects 
for a future detection of this parameter. In Section II 
we review the formalism needed to describe perturbation 
evolution of a general fluid with equation of state and 
sound speed. In Section III we review the main observa- 
tional probes of dark energy, and in section IV we discuss 
the current constraints on w and c s . In Section V we dis- 
cuss possible future constraints on the sound speed from 
observations of the CMB, large scale structure, and weak 
lensing. Finally, section VI contains a discussion of the 
results. 
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II. FORMALISM 

The standard treatment of cosmological perturbation 
theory is the now classical paper by Ma and Bertschinger 
|47|. which described perturbation equations for gen- 
eral fluids in both synchronous and conformal gauge. 
The most commonly used numerical Boltzmann solver is 
CMBFAST [42j which uses synchronous gauge through- 
out. For that reason we follow the notation of but 
only derive the perturbation equations in synchronous 
gauge because we use a modified version of CMBFAST 
for all numerical computations. In this gauge the pertur- 
bation equations for a general, imperfect fluid have been 
derived several times in the literature (see for instance 
[HJ E0- EH ) • Here we simply review the equations as 
derived in Ref. [H^ . 

The adiabatic speed of sound for any fluid is given by 



P 



w 



c„ = — = w - 

a p 3H(l + w) 



(1) 



However, for an imperfect fluid the speed of sound is 
given by the more general relation 



SP 
Tp 



(2) 



If c s ^ c a , entropy perturbations will develop, deter- 
mined by 



P \ P 



(3) 



While many components in themselves behave as perfect 
fluids (dust, radiation, etc), the combination of two per- 
fect fluids with different speeds of sound can produce an 
imperfect fluid in which c s ^ c a . 

Many models of dark energy can also be described as 
imperfect fluids. For example a single scalar field model 
always has c s = 1, while c can be very different. 

In synchronous gauge the density and velocity pertur- 
bations take a particularly simple form [5^| 



-(1 + w) 



k 2 



9 = -H(l-3c 2 s 



1 + w 



5 - k 2 a, 



(4) 



where a is the anisotropic stress. This term arises from 
the next order in the Boltzmann hierarchy. In most cases 
this term can be neglected, and the Boltzmann equation 
can be truncated to form a closed system of equations, 
similar to how the Enskog expansion is performed |5.l| . 
These equations are then straightforward to solve, pro- 
vided that we know w and c 2 s . It should be noted that 
the above equations are valid even for time- varying w. 

However, our approach will now be to take w to be 
constant, as has been done in many previous analyses 



[5J, |§a, |§a, |5j| • At present there is no indication that w 
is varying. Even though the present Type la supernova 
data seem to favour a rapid evolution of w, this indication 
vanishes if all available cosmological data is analysed 
l6lL |6"2 (for other discussions of a time- varyin g w , see for 
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c 2 is also taken to be a constant in our analysis. Al- 
though it is possible to produce models with extremely 
rapidly varying speed of sound, such as some fc-essence 
models, as it will turn out the present constraints on c 2 
are so weak that adding parameters to describe a time 
variation of c 2 is not justified at present. 

Finally, it should be noted that the analysis can be 
carried to it; < —1 without any problems, as long as w 
is constant and c 2 > 0. In that case there are no sin- 
gularities in the perturbation equations and they can be 
solved without any problems. We therefore also include 
this region in the likelihood analysis in the next section. 

From the perturbation equations it is clear that for w — 
— 1 there can be no evolution of perturbations in the dark 
energy since 6 = at all times (for the possibility of initial 
dark energy perturbations, see for instance |58tl59| ). 

Furthermore, if w is sufficiently negative, dark energy 
only dominates at z ~ 0. This means that the speed of 
sound can only affect the very largest observable scales 
for these models. 



III. OBSERVATIONAL PROBES 
A. The cosmic microwave background 

The temperature fluctuations are conveniently de- 
scribed in terms of the spherical harmonics power 
spectrum Ct,i = (|a;m| 2 ), where ^?-(#, 4>) = 
S; m a im,Yim(0, 4>) ■ Since Thomson scattering polarizes 
light, there are also power spectra coming from the po- 
larization. The polarization can be divided into a curl- 
free (E) and a curl (B) component, yielding four inde- 
pendent power spectra: Ct,i, Ce,i, Cb,i, and the T-E 
cross-correlation Cte,i- 

The speed of sound of dark energy mainly affects the 
CMB spectrum at the very largest scales via the late 
Integrated Sachs Wolfe effect. In Fig. ^ we show CMB 
spectra for several different models. From the figure it is 
evident that there can only be significant changes when 
w is close to 0. It should also be noted that when there 
are no dark energy perturbations (as in the case where 
the acceleration is caused by modified gravity or other 
non-fluid like effects) the CMB spectrum in general is 
different. From the perturbation equations it is clear that 
there is no simple solution where 6 = at all times unless 
w = —I. The reason is that, independent of c 2 , any 
fluid is affected by the metric term h. In principle it will 
therefore be possible to distinguish any model in which 
the dark matter behaves as a fluid from models where 
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FIG. 1: The CMB spectra for three different values of w. The FIG. 2: The matter spectra for three different values of w. 



full lines are for c 2 a — 0, the dashed for ci = 
line is without dark energy perturbations. 



1, and the dotted 



The full lines are for c s = 0, the dashed for c a = 
dotted line is without dark energy perturbations. 
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there are no dark energy perturbations. However as w 
approaches -1 the difference vanishes. 



B. Large scale structure 

Any large scale structure survey measures the corre- 
lation function between galaxies. In the linear regime 
where fluctuations are Gaussian the fluctuations can be 
described by the galaxy-galaxy power spectrum alone, 
P(k) — \Sk, gg \ 2 ■ In general the the galaxy-galaxy power 
spectrum is related to the matter power spectrum via a 
bias parameter, b 2 = P gg /P m . In the linear regime, the 
bias parameter is approximately constant, so up to a nor- 
malization constant P gg does measure the matter power 
spectrum. 

In Fig. [5] we show matter power spectra for several 
different models. Again, the differences are negligible as 
soon as w approaches —1. 



C. Weak lensing 

At present there are no large scale weak lensing sur- 
veys. However, in the future weak lensing will be one of 
the main probes of cosmology. The advantage of mea- 
suring weak lensing is that it measures evolution at late 
times where dark energy contributes significantly to the 



energy density. In a flat universe the angular convergence 
power spectrum of weak lensing is, in the normally used 
Limber approximation, given by |81| 



p, = -fftn 



XH 



9 2 (X) 
a 2 (\ 



Ki ( ~>X ) (l \ 



(5) 



Here, Hq is the present Hubble parameter, fl m is the 
matter density, I is the multipole moment, and x is the 
radial distance. P n \ is the non-linear matter power spec- 
trum at wave number l/x and redshift corresponding to 
a horizon distance of X- g(x) is the lensing probability 
function 



XH 



n(x') 



X ~x 
X' 



dx- 



(6) 



n(x) is the distribution function of sources, normalised 
to J n(z)dz = 1. Since sources exist at a distribution 
of redshifts which are measurable, it is even possible to 
divide them into redshift bins and do dark energy to- 
mography. The CMB itself can also be used as a source 
in which case n can be taken as a delta function with 
z ~ 1100. We have calculated lensing spectra for the 
same models as before, assuming that all sources are lo- 
cated at z = 2. The non- linear power spectrum was 
calculated using the Peacock-Dodds formalism. This ap- 
proach is semi-analytic and may not be accurate enough 
for real parameter extraction. However, for calculating 
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FIG. 3: The lensing spectra for three different values of w. 



The full lines are for c 3 — 0, the dashed for c s = 
dotted line is without dark energy perturbations. 
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fi = fl m +Q X 1 Fixed 

h 0.72 ± 0.08 Gaussian [83] 

n b h 2 0.014-0.040 Top hat 

n a 0.6-1.4 Top hat 

r 0-1 Top hat 

Q Free 

b Free 



TABLE I: The different priors on parameters used in the like- 
lihood analysis. 



Hubble parameter, and r, the optical depth to reioniza- 
tion. The normalization of both CMB and LSS spectra 
are taken to be free and unrelated parameters. The priors 
we use are given m Table D 

Likelihoods are calculated from \ 2 so that in 2- 
dimensional plots the 68% and 95% regions are formally 
defined by A% 2 = 2.30 and 6.17 respectively. Note that 
this means that the 68% and 95% contours are not nec- 
essarily equivalent to the same confidence level for single 
parameter estimates. 



A. Supernova luminosity distances 



differences between models it is appropriate, as discussed 
in Ref. jl3 ■ The results are shown in Fig. |3| As before 
the differences are significant only when w is close to zero. 
However, as will be discussed in section V, future lensing 
surveys will improve the sensitivity to the sound speed 
of dark energy significantly. 



We perform our likelihood analysis using the "gold" 
dataset compiled and described in Riess et al consist- 
ing of 157 SNIae using a modified version of the SNOC 
package j8j|. 

B. Large Scale Structure (LSS). 



D. Type la supernovae 

Measurements of the cosmological relationship be- 
tween luminosity distance and redshift in itself has ho 
sensitivity to the speed of sound of dark energy because 
it only probes the evolution of the homogeneous back- 
ground. However, since these measurements significantly 
constrain other cosmological parameters they break some 
of the degeneracies between c s and other cosmological pa- 
rameters. 



At present there are two large galaxy surveys of com- 
parable size, the Sloan Digital Sky Survey (SDSS) H 
and the 2dFGRS (2 degree Field Galaxy Redshift Sur- 
vey) 0. Once the SDSS is completed in 2005 it will be 
significantly larger and more accurate than the 2dFGRS. 
In the present analysis we use data from SDSS, but the 
results would be almost identical had we used 2dF data 
instead. In the data analysis we use only data points 
on scales larger than k — 0.15ft./Mpc in order to avoid 
problems with non-linearity. 



IV. CURRENT CONSTRAINTS 

Using the presently available precision data we have 
performed a likelihood analysis for the two dark energy 
parameters w and c s . 

As our framework we choose the minimum standard 
model with 6 parameters: f2 m , the matter density, the 
curvature parameter, f^, the baryon density, Hq, the 



C. Cosmic Microwave Background. 

The WMAP experiment has reported data on Ct,i and 
Cte,i as described in Refs. 0, |j, |8y, |83, Hg. We have 
performed our likelihood analysis using the prescription 
given by the WMAP collaboration 0, 11 IM Hi El which 
includes the correlation between different Cj's. Fore- 
ground contamination has already been subtracted from 
their published data. 
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FIG. 4: The 68% (dark) and 95% (light) likelihood contours 
for w and c 2 3 for WMAP and SDSS data. 

D. Results 

In Fig.0]we show the results of thew likelihood analysis 
for WMAP and SDSS data only. As is also found in 
many other analyses of the same data for w alone, the 
data is compatible with uu = — 1. Furthermore there is 
no constraint on the speed of sound. Note here that the 
bound in w comes from the combination of CMB and LSS 
data. If only CMB data is used then very high values of 
w are allowed. This in turn means that as w increases the 
sensitivity to c 2 s also increases. This effect was found in 
the WMAP data by Bean and Dore [52| , and if only the 
WMAP data is considered there is a tentative indication 
that c 2 < 1. However, as soon as LSS data is added the 
allowed high w region disappears and there are significant 
constraints only on w. 

Next, Fig. [5] shows the same analysis, but only using 
SNI-a data. These measurements are only sensitive to w 
so likelihood contours are vertical lines. 

Finally, in Fig. we show the likelihood analysis com- 
bining all current constraints. As can be seen there is 
again no sensitivity to c 2 because the constraint on w is 
so strong that it rules out the region where changes in 
the speed of sound cause significant effects. 

1. Comparison with models with no dark energy 
perturbations 

However, it is also interesting to compare the likelihood 
analysis where the dark energy is a fluid with well-defined 
equation of state and speed of sound to one where there 
are no perturbations in the dark energy. As mentioned 
before there if no continuous transition between models 
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FIG. 5: The 68% (dark) and 95% (light) likelihood contours 
for w and c\ for SNI-a data. 
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FIG. 6: The 68% (dark) and 95% (light) likelihood contours 
for w and c 2 s for WMAP, SDSS, and SNI-a data. The vertical 
full lines are 68% contours for the model with no dark energy 
perturbations, and the vertical dashed lines are 95% contours 
for the same model. 



with any finite c s and models with no dark energy per- 
turbations because any fluid is subject to evolution in S 
from the metric term h in the perturbation equations. 
The difference between models with no dark energy per- 
turbations and fluid models is therefore also quite inter- 
esting to probe. The vertical lines in Fig. [fj] show the 
likelihood analysis with the same data, but without dark 
energy perturbations. As can be seen the allowed region 
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is shifted to slightly higher values of w. The best fit 
X 2 is 1626.1 for the model with no dark energy pertur- 
bations, whereas it is 1625.5 for the fluid model. Con- 
sidering that the fluid model contains one more fitting 
parameter, there is no evidence at present for or against 
perturbations in the dark energy. The model with no 
perturbations has 1516 degrees of freedom, whereas the 
fluid model has 1515. This means that x 2 /d.o.f = 1.073 
in both cases. 



where = (CI, ftj,, H , n s , r, . . .) is a vector describing 
the given point in model parameter space and C(9) = 
(xx T ) is the data covariance matrix. In the following we 
shall always work with data in the form of a set of power 
spectrum coefficients, Xi, which can be either C'i or P(k). 

If the data points are uncorrelated so that the data 
covariance matrix is diagonal, the likelihood function can 
be reduced to C oc e~ x 1 , where 



V. FUTURE MEASUREMENTS 

While there are no current constraints on the dark 
energy speed of sound, measurements will improve sig- 
nificantly in the future. We have performed a simple 
Fisher matrix error estimation of the performance of fu- 
ture measurements. The estimate is along the lines of 
that presented in Ref . , but includes future large scale 
structure surveys. 

It is possible to estimate the precision with which the 
cosmological model parameters can be extracted from a 
given hypothetical data set. The starting point for any 
parameter extraction is the vector of data points, x. This 
can be in the form of the raw data, or in compressed form, 
typically the power spectrum (C; for CMB and P(k) for 
LSS). 

Each data point has contributions from both signal and 
noise, x = x s i gna i -I- x no i SC . If both signal and noise are 
Gaussian distributed it is possible to build a likelihood 
function from the measured data which has the following 
form l9Hl 



1 



£(6)cxexp - ^[C(e)- 1 ] 



(7) 



(^i,obs theory) 



i=l 



a(Xi) 2 



(8) 



is a x 2 -statistics and iV max is the number of power spec- 
trum data points [9l| . 

The maximum likelihood is an unbiased estimator, 
which means that 



<e> = e„. 



(9) 



Here 0o indicates the true parameter vector of the under- 
lying cosmological model and (9) is the average estimate 
of parameters from maximizing the likelihood function. 

The likelihood function should thus peak at ~ Oo, 
and we can expand it to second order around this value. 
The first order derivatives are zero, and the expression is 
thus 



\ 2 ■ 

A mm 



1,3 



dx k dx k 



,obs *£/c) 



d 2 x k 



(9j - 9), 



(10) 



where i,j indicate elements in the parameter vector O. 
The second term in the second derivative can be expected 
to be very small because (xfc,obs — %k) is in essence just 
a random measurement error which should average out. 
The remaining term is usually referred to as the Fisher 
information matrix 



inequality 90] that 



n{9 l ) = y/(F~*Ti 



(12) 



for an optimal unbiased estimator, such as the maximum 
likelihood lal. 



F;, 



d 2 x 2 
89,89, 



1 8x k 8x k 



S a(x k )2~89~d9~ 



(11) 



The Fisher matrix is closely related to the precision with 
which the parameters, 0j, can be determined. If all free 
parameters are to be determined from the data alone 
without any priors then it follows from the Cramer- Rao 



A. Mock LSS surveys 

For purposes of parameter estimation the most impor- 
tant parameter in galaxy surveys is the effective volume, 
defined as 



Veff = 



n(r)P(k) 



l+h(r)P(k) 



rf 3 / 



(13) 
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In the above equation n(r) is the selection function. The 
simple interpretation of V is that it is the volume avail- 
able for measuring power at wavenumber k. 

In the following we shall assume that the survey is vol- 
ume limited, meaning that the selection function is con- 
stant throughout the survey volume. If the survey is flux 
limited the selection function is much more complicated. 
In the region where P(k) > 1/n, V e g is independent of k 
and equal to the total survey volume. 

Essentially this means that, with certain restrictions, 
we can use just one free parameter, V c g, to describe a 
hypothetical galaxy survey. 

It was shown in Ref. [93j that the contribution to the 
Fisher matrix from such a galaxy survey can be written 
as 



Ft, 



2tt 



dlnP(fc) dlnP(fc) 



89, 



89, 



w(k)dlixk, (14) 



where the weight-function is w(k) = V c g/X 3 and A = 
2ir/k. The upper limit of the integral should be taken to 
be A; cu t, which is the scale where non-linearity starts to be 
significant. In principle the lower limit, fc m in, should be 
zero but at large scales the assumption that P(k) ^ 1/n 
breaks down. However, by far the most of the weight in 
the above integral comes from k close to the upper limit. 
Therefore, as long as the k where P(k) = 1/n is much 
smaller than fc max the error from taking V e s and k min = 
is quite small. 

It should be noted that the above integral expression 
is quite crude. However, it offers a very simple way of 
estimating parameter estimation errors from galaxy sur- 
veys. The error arising from it can be of order a factor 2, 
leading to an error in the estimated a(9i) of order 2 1 / 2 . 

Instead of V e g we use A e ff = (3T4ff/47r) 1 / 3 as the free 
parameter. As discussed in Ref. [93] the SDSS Bright Red 
Galaxy (BRG) survey has an effective volume of roughly 
(l/i -1 Gpc) 3 , corresponding to X e s — 620/i _1 Mpc. 

Note that the number of independent Fourier modes 
on a given scale, k, enclosed within the survey volume is 
proportional to V e g . Therefore it essentially corresponds 
to the factor (21 + 1) for the CMB measurements which 
measures the number of m- modes for a given In that 
sense both V e g and (21 + 1) are a measure of the lack of 
ergodicity in the given data set. 



B. Mock CMB experiments 

For a mock CMB experiment we assume it to be cosmic 
variance (as opposed to foreground) limited up to some 
maximum Z-value. This value can, however, be different 
for temperature and polarization detection. Therefore a 
given hypothetical experiment can be described by only 
two free parameters, ir, ma x and Zp imax . For all experi- 
ments it will be the case that ?T,max > lp tmax . 

In this picture the MAP data will be well described 
by ^T.max — 1000 and lp max = 0, and the Planck data 



by !r ; max — 2500 and lp !max — 1500. In some sense 
Planck can therefore be regarded as the "ultimate" CMB 
experiment because is measures all of the power spectrum 
parameter space not dominated by foregrounds. 

The contribution to the Fisher matrix from such a 
CMB experiment is then 



^ dCi.x „ -\ tr , n \dQ,Y 
F v = E 1^^q- Cov ( C l,x,Ci, Y )- 



1=2 XY 

/T.max 



89, 



E 



■Cov {Cit,Ci iT ) ' ,( 15 ) 



i=i i 



89, 



89, 



where X,Y = T, E, TE. 



The covariance matrices are given by [9J 
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1) 


2 




(21 + 
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r 2 



r 2 



r 2 



CltCite 



CleCite 



(16) 
(17) 

1(18) 
(19) 
(20) 
(21) 



It should be noted here that this approximation relies 
on the assumption of 4-7T sky coverage and no pixel noise 
up to the maximum /. Even though these assumptions 
are not realised in any real experiment they are suffi- 
ciently accurate for estimating the parameter estimation 
accuracy of a given experiment. 

Finally, a word of caution on the treatment used here. 
In the Fisher matrix analysis for CMB no window func- 
tion is used, i.e. it is assumed that all Ci values are un- 
corrected. While this is true for I much smaller than the 
resolution limit for a given experiment, it breaks down 
close Z max . Furthermore, since most of the information 
on cosmological parameters, including the neutrino mass, 
comes from high I this is a real concern. However, by 
choosing / max conservatively in the Fisher matrix analy- 
sis this problem of overestimating precision at I close to 
the resolution limit can be minimized. 



C. Weak lensing measurements 

The uncertainty in the angular convergence spectrum 
for a weak lensing survey is given by 



(21 + l)/ sky 



1/2 



P K 



(lit) 



(22) 
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Data 




CMB only 


6.26 


CMB + LSS 620 


5.67 


CMB + LSSaooo 


5.66 


CMB + LSS2000- 


hWL 3.03 



TABLE II: la uncertainty on c 2 s for different cosmological 
data sets. 

where / s k y is the fraction of sky measured. (7^) — 0.4 is 
the intrinsic ellipticity of galaxies and n ~ 6.6 x 10 8 sr" 1 
is the angular density of galaxies. At small scales non- 
linearity invalidates the assumption of Gaussianity inher- 
ent in the Fisher matrix formalism, and for this reason 
we do not use data beyond I = 3000. For a lensing survey 
of this type the only adjustable parameter is / s k y which 
can at maximum be about 0.65. 



D. Results 

In order to calculate estimated la errors on the vari- 
ous cosmological parameters we need to apply the Fisher 
matrix analysis to a specific cosmological model. 

We choose as the reference model the generic ACDM 
model with the following free parameters: f2 m , the mat- 
ter density, the baryon density, Hq, the Hubble pa- 
rameter, n s , the spectral index of the primordial pertur- 
bation spectrum, r, the optical depth to reionization, Q, 
the spectrum normalization, b, the bias parameter, as 
well as w and (? s . The reference model has the follow- 
ing parameters: £l m = 0.3, Ode = 0.7, Qbh 2 = 0.02, 
H = 70 kms _1 Mpc _x , n s = 1, r = 0, w = -0.8, and 
cl = l. 

As the mock data sets we take the following: For CMB 
we assume a cosmic variance limited experiment up to 
I = 2500 for temperature and 1500 for polarization. For 
LSS we use two cases, one with A e ff — 620/i _1 Mpc and 
one with A c ff ~ 2000/i _1 Mpc. In both cases we take 
fccut = 0.1 h/Mpc. For weak lensing we take / s k y = 0.65. 

The estimated la error bars for the various cases are 
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shown in Table ITT1 Even in the best case scenario a{c 2 s ) ^> 
1 which means that there is little hope of measuring the 
actual value of the speed of sound. 

However, it may still be possible to discern models with 
no dark energy perturbations from models with fluid-like 
dark energy. 



VI. DISCUSSION 

We have calculated constraints on the main dark en- 
ergy parameters, the equation of state, w, and the speed 
of sound, c 2 , from current observational data. 

It was found that there is at present no useful con- 
straint on the speed of sound of such a general fluid com- 
ponent. The reason is that current observations constrain 
w to be close to — 1 where the sound speed of dark energy 
has a negligible influence on perturbation evolution. 

We also compared models in which dark energy be- 
haves as a fluid to models with no dark energy perturba- 
tions. Here, it was found that the allowed region for w is 
shifted to slightly higher values. The x 2 /d.o.f, however, 
is equal to 1.073 in both cases, meaning that at present 
there is no evidence for or against perturbations in the 
dark energy. 

With regards to future observations we have performed 
a Fisher matrix analysis of CMB, LSS, and weak lensing 
measurements. Even in the most optimistic case where 
full sky measurements of CMB polarization and weak 
lensing can be performed the projected sensitivity to the 
speed of sound in fluid models is quite poor, and there is 
little hope of measuring it. However, it might be possible 
to discern between models where there are no dark en- 
ergy perturbations and models in which the dark energy 
behaves as a fluid. 

Finally we note that even if the speed of sound can- 
not be measured, it can still bias measurements of other 
parameters. For example, for the present cosmological 
data, treating the dark energy as a fluid leads to a slightly 
different estimate for the dark energy equation of state 
compared to the case where no dark energy perturbations 
are assumed. 
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